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First, we establish which measures of large-scale perturbations are least afflicted by gauge arti- 
facts and directly map the apparent evolution of inhomogeneities to local interactions of cosmological 
species. Considering nonlinear and linear perturbations of phase-space distribution, radiation in- 
tensity and arbitrary species' density, we require that: (i) the dynamics of perturbations defined 
by these measures is determined by observables within the local Hubble volume; (ii) the measures 
are practically applicable on microscopic scales and in an unperturbed geometry retain their micro- 
scopic meaning on all scales. We prove that all measures of linear overdensity that satisfy (i) and (ii) 
coincide in the superhorizon limit. Their dynamical equations are simpler than the traditional ones, 
have a nonsingular superhorizon limit and explicit Cauchy form. Then we show that, contrary to 
the popular view, the perturbations of the cosmic microwave background (CMB) in the radiation 
era are not resonantly boosted self-gravitationally during horizon entry. (Consequently, the CMB 
signatures of uncoupled species which may be abundant in the radiation era, e.g. neutrinos or early 
quintessence, are mild; albeit non-degenerate and robust to cosmic variance.) On the other hand, 
dark matter perturbations in the matter era gravitationally suppress large-angle CMB anisotropy by 
an order of magnitude stronger than presently believed. If cold dark matter were the only dominant 
component then, for adiabatic perturbations, the CMB temperature power spectrum Ce would be 
suppressed 25-fold. 
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I. INTRODUCTION 

Cosmological observations can be valuable probes of 
matter species whose interactions are too feeble to be 
studied by more traditional particle physics techniques. 
Whenever such "dark" species, including dark energy, 
dark matter, or neutrinos, constitute a non-negligible en- 
ergy fraction of the universe, these species leave gravita- 
tional imprints on the distribution of more readily observ- 
able matter, notably, the cosmic microwave background 
(CMB) and baryonic astrophysical objects. In the frame- 
work of a perturbed cosmological expansion, the dark 
species influence the visible matter by affecting both the 
expansion rate, controlled by their average energy, and 
the metric inhomogeneities, sourced by species' pertur- 
bations. 
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The gravitational signatures of the perturbations are 
particularly informative as they are sensitive to the in- 
ternal (kinetic) properties of the dark species. These sig- 
natures are generally dissimilar for species with identical 
background energy and pressure but different dynamics 
of perturbations, e.g., for self-interacting particles, un- 
coupled particles, or classical fields [1-5]. Perturbations 
of dark species can affect the observed CMB and matter 
power spectra by an order of magnitude. For example, 
the gravitational impact of dark matter perturbations 
suppresses the CMB power at low multipoles I < 200 by 
up to a factor of 25 (Sec. VC). 

Despite the value and prominence of the gravitational 
impact of perturbations, most of the existing descriptions 
of cosmological evolution are misleading in matching dy- 
namics of perturbations in dark sectors to features in 
observable distributions. One long-investigated source of 
ambiguities and erroneous conclusions is the dependence 
of the apparent evolution of perturbations and the in- 
duced by them gravitational fields on the choice of space- 
time coordinates (metric gauge), e.g. [6, 7]. Of course, 
the CMB and large scale structure observables should 
be identical in any gauge. However, in various gauges 
the features of the observable distributions may appear 
to be generated in different cosmological epochs and by 
different mechanisms. 

For example, given the standard adiabatic initial con- 
ditions, the perturbation of the CMB temperature 6T/T 
grows monotonically on supcrhorizon scales in the his- 
torically important and still popular synchronous gauge, 
e.g. [8]; ST/T generally remains constant beyond the 
horizon in a single cosmological epoch but changes dur- 
ing the transition to another epoch in the calculationally 
convenient and intuitive Newtonian gauge, e.g. [9, 10]; 
or, 8T/T is strictly frozen beyond the horizon but evolves 
differently since horizon entry in the spatially flat gauge. 
These descriptions naively suggest different separation 
of the presently measured CMB temperature anisotropy 
into the inhomogeneities of primordial (inflationary) ori- 
gin and those generated by the gravitational impact of 
perturbations in various species (CMB, neutrinos, dark 
matter, dark energy, etc.) 

The use of gauge-invariant perturbation variables 1 [6, 
7, 11, 12] does not remove this ambiguity [13]. In- 
deed, the perturbations of species' density and veloc- 
ity or the metric in any fixed gauge can be written as 
gauge- invariant expressions, which may be called "gauge- 
invariant" (though, clearly, non-unique) definitions of 



these perturbations [13]. Hence, the variety of gauge- 
invariant perturbation variables is at least as large as the 
variety of gauge-fixing methods. 

Nevertheless, without contradicting general covari- 
ance, much of this descriptional ambiguity can be 
avoided. We show that in cosmological applications the 
perturbed evolution can be described by variables whose 
change is necessarily induced by a local physical cause. 
Moreover, when the equations of perturbation dynamics 
are presented in terms of such variables, the structure of 
the equations simplifies considerably. The simpler equa- 
tions allow tractable analytical analysis of perturbation 
evolution for realistic cosmological models. 

We impose two requirements on a measure of dynam- 
ical cosmological perturbations: 

I. The dynamics of the measure is determined com- 
pletely by locally identifiable physical phenomena 
within the local Hubble volume; 

II. The measure is universally applicable to all scales: 
from superhorizon (governed by general relativity) 
to subhorizon (governed by microscopic kinetics). 

The goals of this paper are, first (Sees. II and III), to 
provide a full dynamical description of cosmological in- 
homogeneities in terms of measures which satisfy these 
two requirements. Second (Sees. IV and V), to establish 
the origin of features in the CMB angular spectrum, as 
revealed by this more direct description, and to explore 
the utility of these features for probing the dark sectors. 

Most of the best-known formalisms for the dynamics 
of CMB and matter inhomogeneities are formulated in 
terms of perturbations of physical observables, such as 
radiation temperature T(x M ,n) or proper energy den- 
sity p(x M ) [10, 14-17]. When these perturbations are 
evaluated in the Newtonian gauge, nonsingular on small 
scales, then the small-scale dynamics in these formalisms 
does reduce to special-relativistic kinetics and Newtonian 
gravity. Then, however, the perturbation evolution be- 
yond the Hubble scale is nonlocal 2 and subject to gauge 
artifacts. The artifacts are caused by the nonlocality of 
defining the motion of coordinate observers, hence, infer- 
ring the components of tensors and the splitting of dy- 
namical variables into their background values and per- 
turbations from the gauge conditions. While some gauge 



2 For illustration, a perturbation of energy density p in a dust (zero 
pressure) universe evolves in the Newtonian gauge (26) as 



1 A "gauge-invariant perturbation variable" generally cannot be 
measured by a local observer. (For example, the local values of 
the gauge-invariant Bardeen potentials <i> and \f [6] are deter- 
mined by spacctime curvature outside the observer's past light 
cone.) Thus the gauge-invariant perturbations should be dis- 
tinguished from gauge-invariant local observables or from local 
physical tensor quantities, such as the energy-momentum ten- 
sor T^ v or the curvature tensor a . 



5p + 3- 5p + pVi^ 1 
a 



3p<S>, 



with 



<J> = AirGa 2 — 



8p - 3- p— - \7 t v' 
a 



Although the first equation above is obtained by linearization 
of a causal relation T°£ = 0, the operators 1/V 2 in the second 
equation signal that after elimination of <t> the evolution of Sp 
becomes nonlocal. This does not violate causality because <5p 
cannot be measured locally. 
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conditions, e.g. synchronous, can be imposed locally, the 
resulting descriptions tend to be singular, contain gauge 
modes, and fail to reduce to Newtonian gravity on small 
scales. 

Studies of the connection of the observed cosmological 
inhomogcneitics to quantum fluctuations of an irmaton 
field during inflation have led to an extensive list of vari- 
ables which under certain conditions freeze (become con- 
stant in time) beyond the horizon. The best known exam- 
ples are the Bardeen curvature C [6, 18] (conventionally 
interpreted either as a perturbation of intrinsic curva- 
ture on uniform-density hypersurfaccs or as density per- 
turbation on spatially flat hypersurfaces) and the curva- 
ture perturbation 1Z on comoving hypersurfaccs, e.g. [19]. 
Both £ and 1Z are frozen for adiabatic supcrhorizon per- 
turbations. For the general, nonadiabatic, superhorizon 
perturbations, the curvature perturbations £ a [20] on the 
hypersurfaces of uniform energy density of an individ- 
ual minimally coupled perfect fluid a were shown to be 
also frozen ("conserved"). 3 Finally, perturbed cosmolog- 
ical evolution has been described in terms of conserved 
spatial gradients of various quantities [22-24] . The evo- 
lution of any of these variables is more robust to gauge 
artifacts on supcrhorizon scales than, for example, that of 
5T/T, Sp/p, or (the Bardeen) potentials $ and ^ in the 
Newtonian gauge. However, neither the uniform-density, 
nor comoving, nor spatially flat hypersurfaces reduce to 
Newtonian hypersurfaces on subhorizon scales. Nor does 
the evolution of gradients tend to the familiar picture of 
particles and fields propagating in Minkowski spacetime. 

In Sec. II we consider natural measures for nonlinear 
perturbations of phase-space distribution and radiation 
intensity and for linear perturbations of species' density 
that conform to both requirements I and II. Moreover, we 
prove in Sec. II D that, although there is no "physically 
preferred" description of perturbed evolution, the over- 
all change of density perturbations of any species during 
horizon entry is the same in any description which satis- 
fies certain conditions formalizing requirements I and II. 
This change is different in most of the traditional for- 
malisms, violating these requirements. 



the perturbed cosmological dynamics. The proposed ap- 
proach has broad scope of applicability, providing an eco- 
nomical and physically adequate description of phenom- 
ena which involve inhomogeneous evolution of multiple 
species and non-negligible general-relativistic effects. Ex- 
amples include the physics of inflation, reheating, the 
CMB, and cosmic structure. In this paper we focus on 
the last two topics. 

After a concise review of the evolution of perturba- 
tions on superhorizon scales and during horizon entry in 
Sec. IV, we apply in Sec. V the developed formalism to 
study the gravitational signatures of various species in 
CMB temperature anisotropy and large scale structure. 

Contrary to a popular view, e.g. [17, 25-29], based on 
the study of traditional proper perturbations in the New- 
tonian gauge, we find that CMB perturbations in the 
radiation era are not resonantly boosted by their self- 
gravity (Sec. VB 1). As a consequence, the gravitational 
signatures of dark species in the radiation era, such as 
neutrinos or a dynamical scalar field (quintessence [30, 
31]), are moderate; as these species, even when they are 
abundant, do not untune any physical resonant amplifi- 
cation. Fortunately, due to low cosmic variance on the 
corresponding scales and the existence of characteristic 
nondcgcncratc signatures (Sec. VB2), we can still ex- 
pect meaningful robust constraints on the nature of the 
dark radiation sector. 

On the other hand, the gravitational impact of dark 
sectors' perturbations on the CMB in the matter era is 
found to be an order of magnitude stronger than the 
traditional interpretations suggest (Sec. VC). 

We summarize our results in Sec. VI. Appendix A 
presents linear dynamical equations for scalar perturba- 
tions of typical cosmological species in terms of the sug- 
gested measures. Appendix B summarizes the formulas 
for the scalar transfer functions and angular power spec- 
tra of CMB temperature and polarization. The main 
notations to be used in this paper are listed in Table I. 



II. MEASURES OF PERTURBATIONS 
A. Phase-space distribution 



In Sec. Ill we describe a closed formulation of lin- 
ear dynamics of scalar perturbations in terms of the 
above measures evaluated in the Newtonian gauge. We 
find that this formulation has several technical advan- 
tages, ultimately related to its tighter connection be- 
tween the perturbation measures and the causality of 



3 The perturbations f a [20] are conserved for uncoupled perfect flu- 
ids [20] or species whose perturbations are internally adiabatic [4] 
only if gravitational decays [21] of species into another type of 
species arc negligible. This should be a reasonable assumption 
for realistic applications to the evolution of the CMB and large 
scale structure. 



Our first example of a quantity whose dynamics can be 
fully determined by physics in the local Hubble volume 
is the one-particle phase-space distribution of classical, 
possibly relativistic, point particles [photons, neutrinos, 
or cold dark matter (CDM)] f{T,x l ,P t ) = dN/d 3 x l d 3 Pi. 
(By default, Latin indices range from 1 to 3 and Greek 
from to 3.) We consider / as a function of the co- 
ordinate time t = x° , spatial coordinates x % , and, cru- 
cially for the following results, canonically conjugate mo- 
menta Pi. The distribution / evolves according to the 
Boltzmann equation 



/ + 



dx i df 
dr dx i 



mdf_ 

dr dPi 



C 



(1) 
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Symbol 


Meaning 






Definition 


r 


Coordinate time, x u [Conformal time in the FRW backgroun 


d] 




Sec. II A [eq. (26)] 


X 


Spatial coordinates, x 1 [Comoving coordinates in the FRW backg 


round] 


[eq. (26)] 


Pi 


Canonical momenta 






Sec. II A 


P 


(ELi P ') 1/2 in an y metric 






eq. (3) 


rii 


Direction of propagation, Pi/P 






eq. (3) 


fix", Pi) 


Phase-space distribution 






Sec. II A 


df(x\P) 


Canonical perturbation of /, f(x l , Pi) — f(P) 






Sec. II A, eq. (A3) 




Conformal intensity of radiation 






eq. (7) 


b{x^,n t ) 


Perturbation of radiation intensity, I /I — 1 






eq. (9) 


/ck/3, I(X0, Lap 


Describe polarized photons 






Sec. A 3 c 


p(x»), p(x") 


Energy density and pressure of species 






p. 6 


d(x^) 


Perturbation of species' coordinate number density 






eq. (20) or (27) 


(for a fluid of particles, Sn coo /n coo with n coo = dN/d 3 x) 






(P- 7) 




Normal bulk velocity of species 






eq. (17) 


v'^x") 


Coordinate bulk velocity 






eq. (21) 


u(x*) 


Velocity potential of scalar perturbations, Vi = — ViW 






eq. (30) 


o-(V) 


Scalar potential of anisotropic stress 






eq. (32) 




Scalar multipole potentials of i (in particular, do = d, di — i 


, d 2 


= ff ) 


eq. (A8) 




Scalar multipole potentials of photon polarization 






eq. (A21) 




A classical scalar held (quintessence) 






Sec. A 4 




Metric tensor, with the signature ( — , +, +, +) 








a 


Scale factor in the FRW background 






eq. (14) 




General perturbation of the metric, Sg^/a 2 






eq. (14) 


D and e 


General scalar perturbations of the spatial metric gij 






eq. (29) 


$ and $ 


Scalar perturbations of the metric in the Newtonian gauge 






eq. (26) 




Reduced curvature, D + ±V 2 e = -\ a 2 (V" 2 ) (3 lR, of 3-slices 


Pa = 


= const 


Sec. Ill A 


overdot, ' 


d/dr 








H 


Coordinate expansion rate [conformal, in the FRW background] 




Sec. II A [a /a] 


7 


A-kGcl (p + p) in the FRW background 






eq. (34) 


t cf f and 9 ff 


Effective intensity and temperature perturbations 






eqs. (47) and (48) 


Cs 


Speed of sound in the photon-baryon plasma 






eq. (A10) 


R b 


Baryon to photon enthalpy ratio, 3pf,/(4 / o 7 ) 






Sec. V A 


S 


Acoustic horizon, Q c s dr (when Rb <C 1, r/y/i) 






Sec. V B 2 


<P 


Phase of acoustic oscillations, kS(r) 






eq. (50) 




Mean r of a photon collisionless flight 






eq. (A14) 


&T/T\ in 


Primordial superhorizon perturbation of CMB temperature, 




\k<n 


eq. (57) 


AT/T 


Presently observed perturbation of CMB temperature 




eq. (49) 



TABLE I: Summary of the main notations. 



than the temporal scale of local cosmological expansion, 
Hr 1 , where for nonlinear theory Ti. = ln(det g^). 
Then from cq. (2), where both the second and third 
terms contain spatial gradients (df/dx l and g^v.i), we 
see that /// = 0(A _1 ) <C H. In the limit of super- 
horizon inhomogeneities (XH — ► oo) the phase-space dis- 
tribution /(t, x l , Pi) of minimally coupled particles be- 
comes time- independent (frozen). Then its background 
value, /, and perturbation, df = f(x fi ,Pi) — f(P), are 
also frozen. 

Conversely, an apparent temporal change of df (at 
fixed x % and Pi in any regular gauge) can always be 
attributed cither to the gravity of physical inhomo- 
geneities within the local Hubble volume or to local non- 
gravitational interaction. 

The majority of contemporary formalisms for cosmo- 
logical evolution in phase-space, e.g. [10, 15, 17, 33, 34], 
work not with the canonical momentum of particles Pi 



(for a systematic formulation of kinetic theory in gen- 
eral relativity see, e.g., [32].) If the particles interact 
only gravitationally, their canonical momenta P, coin- 
cide with the spatial covariant components of the parti- 
cle 4-momenta P M : P} = gi^P 11 , where g^ v is the metric 
tensor. Then dx l /dr = P l /P a , dPi/dr is given by the 
geodesic equation, and C, accounting for two-particle col- 
lisions, vanishes. Hence, 

J + P o g x i + 2 P0 dPi ~ [ ' 

We stress that this equation and the following observa- 
tion apply to the fully nonlinear general-relativistic dy- 
namics. 

Let us consider the minimally coupled particles in 
the perturbed spacetime, possibly populated by addi- 
tional species. Let us also suppose that in certain co- 
ordinates the spatial scale, A, of inhomogeneities in the 
particle distribution and in the metric is much larger 
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but rather with proper momentum p, measured by coor- 
dinate or by normal observers. These formalisms typi- 
cally consider the phase-space distribution as a function 
of q = ap, the proper momentum rescaled by the back- 
ground scale factor. The evolution of the correspond- 
ing perturbation Sf(x^,q) = f(x^,q) — f(q), as well as 
of perturbations of integrated proper densities and in- 
tensities, depends on the nonlocal procedure for defining 
the observer's frame from the gauge condition. In typ- 
ical fixed gauges 8f{x^,q), unlike df{x^ 1 P i ) 1 generally 
changes even in the absence of local physical inhomo- 
geneity or local non-gravitational coupling. 4 

Thus the stated in the introduction criteria I and II 
for a measure of cosmological inhomogeneities are ful- 
filled for a perturbation of one-particle phase-space dis- 
tribution f{x^ L 1 Pi), provided Pi are the particle canonical 
momenta. A description in terms of phase-space distribu- 
tions is, however, unnecessarily detailed for most cosmo- 
logical applications. Instead, it is more convenient to use 
integrated characteristics of the species, such as inten- 
sities (for radiation of photons or other ultrarelativistic 
species) or energy densities and momentum-averaged ve- 
locities (for arbitrary species, including non-rclativistic 
particles, fluids, or classical fields). In the next two sub- 
sections we discuss the general-relativistic measures of 
perturbations in intensity and density that also satisfy 
criteria I and II. 



B. Radiation intensity 

We start by considering inhomogeneities in the inten- 
sity of any type of cosmological radiation, such as CMB 
photons or relic neutrinos at the redshifts at which the 
kinetic energy of the particles dominates their mass. We 
describe the direction of particle propagation by 

P 3 
m=^ where P 2 = ^/f, (3) 

i=i 

so that 

i>?=i. ( 4 ) 

i=l 

We also introduce no = P^/P (for the ultrarelativistic 
particles g^n^ = 0) and n^ = g» v n v = P^/P. The 



motivation for the noncovariant definitions (3) is to de- 
scribe radiation transport in terms of variables which are 
fully specified by the "dynamical" quantity f(T,x l ,Pi) 
but are unmixed with the gauge-dependent and, in gen- 
eral, noncausally evolving metric g^. 

The energy-momentum tensor of radiation equals 

Substituting P^ — n^P and P^ = n^P, we can rewrite 
it as 

J V=<7 n 

-a directional average of an expression which depends 
only on the local value of the metric and the quantity 

/>OC 

l(x^m)= / P 3 dPf(x», ni P). (7) 

JO 

The variable (7), to be called the conformal intensity 
of radiation, inherits such useful properties of f(x^,Pi) 
as timc-indcpcndcncc for supcrhorizon inhomogeneities 
of non-interacting particles and reduction to the proper 
intensity, dE / (dVcPh), in the Minkowski limit. 

The dynamics of 7(x^,ni) for minimally coupled ra- 
diation in a known metric is given by a closed equa- 
tion [35] which is straightforwardly derived by integrating 
the Boltzmann equation (2) over P 3 dP and remembering 
that for the ultrarelativistic particles g^n^n^ = 0: 5 

This equation confirms that for decoupled particles which 
are perturbed only on supcrhorizon scales (the gradients 
dl/dx 1 and gp, v ,i are negligible) the conformal intensity 
freezes: 1 = 0. 

Now we can describe perturbation of the intensity by 
a variable 

t(a^,ni) = (9) 

where I is a time-independent background value of /. As 
intended, the perturbation i is time-independent in full 
(nonlinear) general relativity for supcrhorizon perturba- 
tions of uncoupled radiation. This variable also reduces 



4 For example, 
gauge (26), g» 



for linear perturbations in the Newtonian 
= Pi/ {I - *) [10, 17, 34]. Hence, 

Sf(x»,q) = df(x»,Pi) -VPi^L. 

aPi 

While the first term on the right-hand side evolves causally on 
supcrhorizon scales, the Newtonian potential "I" does not, and so, 
neither does 8f(x ll ,q). 



5 In the last terms of eqs. (8) and (10), 
with respect to the components of n, = 
condition (4), can be naturally defined as 



d 
drii 



E 



f PdtM a 
V dPi J 9/a 



the partial derivatives 
: Pi/P, constrained by 

a 



where fj, a = (fii,fJ.2) are any two independent variables parame- 
terizing n. 
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to the perturbation of proper intensity in spatially homo- 
geneous and isotropic, or Friedmann-Robertson- Walker 
(FRW), metric. The full dynamical equation for l evolu- 
tion follows trivially from eq. (8): 



2m(i + l) 



1 dt 
2~dn~ 



(10) 



The value of the introduced intensity perturbation t 
is generally gauge-dependent. If spacetime is spatially 
homogeneous and isotropic then the FRW metric is nat- 
urally preferred as the metric manifesting the spacetime 
symmetries. In this metric i is unambiguous and co- 
incides with the perturbation of proper intensity. As 
shown next, the value of i is also unambiguous for su- 
perhorizon perturbations: it is the same in any of the 
gauges which provide a nonsingular description of the 
supcrhorizon evolution. (More general arguments for the 
uniqueness of "physically adequate" mapping of super- 
to subhorizon perturbations are given in Sec. II D.) 

Under an infinitesimal change of coordinates x^ — > 



Sx^ the perturbation i transforms as 



l(S^, hi) = L(x^,ni) - iriin^Sx^, 



where by eq. (3) 



= Ui- (Sij - mrij)n^6x 



(11) 



(12) 



Suppose that in at least some gauges the superhorizon 
evolution of cosmological perturbations appears nonsin- 
gular (regular) . Namely, that in these gauges the pertur- 
bations of all the dynamical and metric variables Q have 
a moderate magnitude \5\ = \SQ/Q\ < 1 and are homoge- 
neous on the spatial scales less than A ~ |<V^,i| 3> H^ 1 . 
A transformation between two gauges either of which is 
regular, as defined here, has to be restricted as 



1^1 < w~ 



and 



1^1 < n-'/x. (is) 



Given these restrictions and the result di/dx^ = 0(A _1 ), 
following from them and from eq. (10), the values of rij 
and l are unchanged by the gauge transformation (11)- 
(12) in the superhorizon limit Hr 1 /A — > 0. 6 

The equations given so far applied to the full general- 
relativistic dynamics. We conclude this subsection by 
presenting their linearized versions. We consider the per- 
turbed FRW metric 



(14) 



Most of the gauge-dependent measures of perturbations vary in 
the H~ 1 /X -* limit even within the restricted class of regu- 
lar gauges. For example, the metric perturbation Sg^ v changes 
under a restricted gauge transformation i M = x^ + Sx 11 , with 
l<5<«l < - 0, by -C Sx g„ v = -(<? M <5° - S°g x „)8x* - 

g^Sx 0. 



with ?7 M „ = diag(— 1, 1, 1, 1). The linear in perturbation 
terms of eq. (10) are 

i + mi t i = 2n i n t *n"h tll ,^, (15) 

where n M = (l,n,). The intensity perturbation i deter- 
mines linear perturbation of radiation coordinate den- 
sity d, introduced in the next subsection, as [eq. (25)] 

d=|(t> ft , (16) 

where ()„ denotes J d 2 ni/4ir. Likewise, i specifies normal 
bulk velocity of the radiation 



P + P' 



where p = -T ° and p = £\ T i/3, as 

Vi = - (njt)ft, 
which follows from eq. (6). 



(17) 



(18) 



C. Energy density 

Next we consider a measure of perturbation of energy 
density for arbitrary species that, again, is directly con- 
nected to local physical dynamics and is universally ap- 
plicable to all scales. While we lift all the restrictions 
on the nature of the species, we limit our construction 
to the linear order of cosmological perturbation theory. 
For the discussion of densities it is sufficient to address 
scalar perturbations and disregard uncoupled vector and 
tensor modes. 

Suppose that the studied, possibly self-interacting, 
species (labeled by a subscript a) do not couple non- 
gravitationally to and are not created in significant 
amount by gravitational decays of the remaining species. 
Such species, to be called "isolated", can be assigned 
a covariantly conserved energy-momentum tensor T£ v . 
In an arbitrary gauge, the covariant energy conservation 
f-aw = in linear perturbation order gives 

PaSpa-paSpa ^ 



(Pa+Pa) 2 



Here we introduced coordinate number density pcrturba 
tion of the species 

Sp a 



d a = 



3D, 



(20) 



Pa + Pa 

where D = i <5(det gtj)/ det g^ = g^/ijj is the spatial 
dilation of the metric (14), and coordinate bulk velocity 
of the species 



T'i 
aO 



Pa +Pa 



(21) 



If the species a are classical particles then v'* gives 
their momentum-averaged velocity measured by coordi- 
nate observers, moving along the lines x l = const. 
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The justification for calling d a "coordinate number 
density perturbation" is the following. Consider iso- 
lated species which are described by an equation of state 
p = p(p) and can be assigned a conserved number N. 
For them, by energy conservation, Sp/(p + p) gives the 
relative perturbation 5n pr /n pr of species' proper num- 
ber density n pr = dN/dV pr . Then d a equals the per- 
turbation 5n coo /n coo of species' coordinate number den- 
sity n coo = dN/d 3 x, where the proper and coordinate in- 
finitesimal volumes are related as dV pl = a 3 (l + 3D)d 3 x. 

The right-hand side of eq. (19) is a gauge-invariant 
quantity, proportional to the so called "nonadiabatic 
pressure" Sp a — (p a / p a )$Pa- It vanishes trivially when- 
ever p a is a unique function of p a - Even if p a and p a are 
not functionally related, for example, if a is a field or a 
mixture of more basic interacting species, the right-hand 
side of eq. (19) vanishes whenever all the proper distribu- 
tions describing the species a are unperturbed in certain 
coordinates. Indeed, then p a Sp a —PaSp a vanishes in those 
coordinates and so, by its gauge invariance, vanishes in 
any coordinates. We will call the superhorizon pertur- 
bations of the species a which appear unperturbed on a 
certain hypcrsurface internally adiabatic. We also note 
that the right-hand side of eq. (19) can be measured by 
a local observer, provided he accepts a background equa- 
tion of state p a = Pa(Pa) which defines the nature of the 
unperturbed species. 

According to eq. (19), the density perturbation d a 
is frozen for superhorizon internally adiabatic perturba- 
tions of isolated species. This superhorizon conservation 
can be easily understood for particles with a conserved 
number, when d a is literally the perturbation of coordi- 
nate density of this number. As long as the chosen gauge 
is nonsingular in the superhorizon limit Hr 1 /A — > 0, 
the particle number n coo in a unit coordinate (but not 
proper!) volume changes infinitesimally over a time inter- 
val At = H-\ by An coo /n coo = At = O^/X). 
Hence, d a = Sn coo /n coo remains frozen. 

Similarly to the perturbation i of radiation inten- 
sity, the coordinate density perturbation d a is generally 
gauge-dependent. Yet, in the superhorizon limit its value 
in the regular gauges is also unambiguous. This is evident 
from the d a transformation 

d a =d a - Sx^ (22) 

under a change of coordinates x^ = x^ + Sx^ , restricted 
for the regular gauges by conditions (13). 

The traditional measures of density perturbation such 
as the perturbations of proper density or 3-curvature are, 
likewise, gauge-dependent. But, unlike d a , those pertur- 
bations vary in the superhorizon limit even among the 
regular gauges. For example, a gauge transformation law 



5 Pa 



5 Pa 



3HSx° 



(23) 



Pa +Pa Pa+ Pa 

shows that proper overdensities Sp a /(p a + p a ) or Sp a /p a 
remain gauge-dependent in the superhorizon limit among 
the regular gauges, restricted by conditions (13). 



As noted in Sec. I, the values of the gauge- variant den- 
sity, curvature, or other perturbations in any fixed gauge 
define certain "gauge- invariant" perturbations [6, 7, 11, 
12], such as Bardeen's e a or £ [6, 18]. The value of d a 
in a fixed gauge also defines a gauge-invariant variable. 
In particular, d a of the Newtonian gauge (Sec. Ill) cor- 
responds to the gauge-invariant expression (28). 

We can simply relate the perturbation of the coordi- 
nate particle number density d with the discussed ear- 
lier perturbation of phase-space distribution df(x^,Pi) 
for classical point particles and with perturbation of in- 
tensity l for ultrarclativistic particles. As follows from 
eqs. (20) and (5), in linear order, 



Jd 3 P l E(P)df 
a 4 (p + p) 



(24) 



where E(P) = (P 2 +m 2 a 2 ) 1 / 2 . For ultrarelativistic par- 
ticles this equation and eqs. (7) and (9) give 



d = - (i) fi . 



(25) 



Since the variable d a is constant in a regular gauge 
for superhorizon internally adiabatic perturbations of iso- 
lated species, any change of d a can be linked to an objec- 
tive dynamical cause within the local Hubble volume. We 
show next that the variable d a remains useful on all scales 
and that the suggested by it relation between density per- 
turbations on super- and subhorizon scales is uniquely 
preferable on physical grounds. 



D. Uniqueness of superhorizon values 

Various measures of scalar perturbations that freeze 
beyond the horizon, for example, the uniform-density or 
comoving curvatures £ [6, 18] or 1Z [6, 7, 19] or the gradi- 
ents of Refs. [23, 24], differ in their evolution on the scales 
comparable to the horizon and smaller. Not all of these 
measures of perturbations remain meaningful in the sub- 
horizon Minkowski limit. Those that do, e.g. the coordi- 
nate density perturbations d a in various regular gauges, 
may still differ from each other on scales comparable to 
the Hubble scale. 

We now prove that, nevertheless, the value of a super- 
horizon density perturbation is the same in terms of any 
measure which 

I'. Freezes for superhorizon internally adiabatic per- 
turbations of uncoupled (isolated) species, 

and additionally, 

II'. Gives the perturbation of species' number density, 
Sn a /n a = 5p a /(pa +Pa), in an Y spatially homoge- 
neous and isotropic (FRW) metric. 

Condition II' is a formal statement of our desire for a 
measure of density perturbation to be such on all scales 
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whenever the concept of density perturbation is unam- 
biguous (c.f. requirement II in Sec. I). We match the 
measure to 6p a /(p a +p a ), rather than, e.g., to 5p a /p a , 
to be consistent with condition I', which then applies au- 
tomatically to an FRW metric. Finally if the dynamics 
of the measure is fully specified by the physics within 
the local Hubble volume (requirement I in Sec. I), then 
the local dynamics of supcrhorizon perturbations should 
be the same whether the metric is globally FRW or not. 
Hence, we impose condition I' for the general metric. 

The proof is applicable to arbitrary isolated species a. 
To be specific, we consider a decelerating universe, where 
perturbations evolve from super- to subhorizon scales. 
We imagine a cosmological scenario in which the initially 
perturbed geometry becomes homogeneous and isotropic 
while the studied scales A are superhorizon (X/H^ 1 ^> 1) 
and remains homogeneous and isotropic ever since. We 
suppose that the chosen gauge provides a regular descrip- 
tion of superhorizon perturbations and leads to the FRW 
metric in the unperturbed geometry. Let d a and d' a be 
two arbitrary measures that both are frozen on super- 
horizon scales and equal to Sp a /(p a + Pa) in the FRW 
metric. Since in our imaginary scenario the metric be- 
comes FRW on superhorizon scales, d a = d' a at all times. 
In particular, d a = d' a beyond the horizon. We now note 
that if the measures describe an instantaneous state of 
the perturbations then the equality d a = d' a retains its 
validity on superhorizon scales regardless of the subse- 
quent perturbation evolution, i.e., regardless of the cho- 
sen cosmological scenario. 

To conclude the proof, we demonstrate that the in- 
voked imaginary scenario is physically realizable, how- 
ever unnatural this realization may be. Given arbitrary 
superhorizon perturbations of the existing cosmological 
species and the metric, we can force them to evolve to a 
homogeneous geometry by adding new minimally coupled 
species, X, which first contribute infinitesimally to the 
overall energy-momentum density. Let these imaginary 
species obey an equation of state vox = Px / Px = — 1/3 
until they become the only dominant component. (The 
choice w = —1/3 corresponds to the borderline between 
deceleration and acceleration, hence, time-independent 
X/H^ 1 .) We set the coordinate number density pertur- 
bation dx of the species X to vanish on superhorizon 
scales, e.g., by our choice of the initial distribution of 
these species. (This requirement is not dynamically nat- 
ural. Yet, due to supcrhorizon conservation of dx, it is 
self-consistent and in principle realizable.) Then, once 
the unperturbed species X dominate, they generate a 
homogeneous and isotropic geometry. The subsequent 
cosmological expansion may be designed to decelerate 
again, e.g., by taking wx at any XjliT 1 equal w at the 
same XjliT 1 in the original model, without the species X. 
Then in the modified model the fictitious species X re- 
main dominant and the geometry unperturbed, as de- 
sired. We stress that this scenario and the species X are 



introduced only as a gedanken experiment. They arc not 
expected to exist in nature. 7 

The above arguments do not imply that the mea- 
sures d a and d' a are identical. For example, the coordi- 
nate number density perturbations (20) differ in different 
gauges, c.f. eq. (22). The identical are their values for su- 
pcrhorizon perturbations and, trivially, for unperturbed 
geometry, when they all give 5p a /(p a + p a )- 

III. DYNAMICS IN THE NEWTONIAN GAUGE 

Any convenient gauge can be used to study cosmolog- 
ical evolution. Nevertheless, of the most popular gauges, 
such as the synchronous, Newtonian, uniform density, 
comoving, or spatially flat, only the Newtonian gauge 
reduces to Newtonian gravity on subhorizon scales. Be- 
cause of this helpful feature, during our further discussion 
of scalar perturbations we impose the Newtonian gauge 
conditions, namely, require zero shift and shear of the 
metric: g^ = and gij oc Sij. For simplicity, we as- 
sume a spatially flat background geometry. Then the 
perturbed metric reads 

ds 2 = a 2 (r) [-(1 + 2$)dr 2 + (1 - 2^)dx 2 } . (26) 

(It is not a priori obvious that the Newtonian gauge pro- 
vides a nonsingular description of superhorizon scales. 
We will return to this subtlety at the end of the section.) 

A. Densities 

In the Newtonian gauge the coordinate number density 
perturbation (20) is 

da = -^--3% (27) 

Pa + Pa 

and the coordinate and normal bulk velocities of species 
coincide: v'* = Vi a [eq. (21)]. 

The Newtonian-gauge variable (27) can be written in 
a gauge-invariant form 

4 Nowt) = ^^+3^ + V 2 e, (28) 

Pa +Pa 

where D and e parameterize the scalar perturbations of 
the spatial part of the metric as 

5g l3 = a 2 [2D <% - 2(V i V i - V 2 ) e}. (29) 



7 Since the appearance of the first astro-ph version of this pa- 
per (astro-ph/0405157) the suggested scenario has been consid- 
ered by several authors [36, 37] for damping scalar perturbations, 
hence, enhancing the ratio of tensor to scalar modes after infla- 
tion. Ref. [38] , however, points out that the required superhori- 
zon perturbations of the species X are highly unnatural to be 
produced dynamically in the real world. 
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Evaluation of the gauge-invariant combination on the 
right-hand side of eq. (28) in other gauges shows that: 

a) c?i Ncwt ^ = 3£ a , where ( a [20] is the (reduced) cur- 
vature Ca = D + ±V 2 e = -\ a 2 (V~ 2 ) (3) i? on the 
hypersurfaces of uniform density of species a (on 
which p a {x) = const); 

b) <ii Ncwt) equals 5n apT /n ap i = Sp a /(p a + Pa) in the 
spatially flat gauge (Sgij = 0); 

c) <ii Ncwt) coincides with d a = Sn acoo /n acoo of 
eq. (20) in any gauge with zero shear (e = 0), re- 
gardless of the second gauge-fixing condition. 

The evolution of density perturbations d a in any gauge, 
including the Newtonian, is described by eq. (19). To 
complete the formalism, we provide the Newtonian gauge 
equations for the evolution of species' velocity, pressure, 
etc. and for the metric perturbations. 



B. Velocities 

For scalar perturbations, we introduce a velocity po- 
tential u a such that 



Vai = -VjU a . 



(30) 



Covariant momentum conservation T^ ;A1 = for isolated 
species gives 

U a = - + $ 



S Pa Aga + 2 v a (31) 
Pa+Pa 



where a scalar potential a a parameterizes the species' 

i — rpi 1 zi r 

ij — 1 aj 3°j" J 



anisotropic stress j = T^j — ^S^T^k as 



Pa + Pa 



(32) 



The pressure perturbation 8p a and anisotropic stress 
potential a a arc determined by the internal dynam- 
ics of the species. In Appendix A we give the com- 
plete Newtonian-gauge equations for the linear evolution 
of scalar perturbations of typical cosmological species: 
CDM, decoupled neutrinos, partially polarized photons, 
baryons, and a classical scalar field. In particular, for 
the following discussion, we note that the perturbation 
of intensity of decoupled radiation (neutrinos or CMB 
photons after their last scattering) evolves in the Newto- 
nian gauge as 



i + mViL = -4n 4 V l ($ + *), 
as follows from eq. (15). 



(33) 



C. Metric 

The Newtonian gravitational potentials $ and ^ do 
not represent additional dynamical degrees of freedom 
but are fully determined by matter perturbations. To 
write the corresponding equations [4], we introduce the 
reduced background enthalpy density 



j = 4irGa 2 (p + p) = 3(i+!£) W 2 



(34) 



(with p = J^Pa, P = J^Pa and w = p/p), species 
enthalpy abundances x a = (p a + p a )/(p + p), and 
enthalpy-averaged perturbations 



5 P 
P + P 



3-D = x a d a , 



(35) 



u = ^x a u a and a = ^x a a a [for the latter two, 
v i = T?/(p+p) = -V i u and (Tj - ±5^) / (p + p) = 
(VjVj — |5yV 2 )cr]. Then the linearized Einstein equa- 
tions [9, 10] become 



i(* - $) = 7(7. 

In real space eq. (36) is solved by 

-r'/ra 



(36) 
(37) 



*(r) = -£ MV 



(d + 3Hu)\ r+r ,, (38) 



with r G = 1/V37 = lij^) 112 ^ 1 . For subhorizon per- 
turbations the term Tiu is negligible and e~ r l TG pa 1. 
Then the induced potential ^(r) obeys the usual 1/r 
Newton's law (and so does However, on large, 

compared to r G , scales ^ is fully determined by the lo- 
cal coordinate overdensities and velocity potentials of the 
matter species: 



Hu 



(for r > r G ), (39) 



as follows from either eq. (38) or eq. (36). The corrections 
to relation (39) from matter inhomogeneities beyond the 
distance r G are suppressed exponentially. 

The presented dynamical equations for the coordi- 
nate measures of density and intensity perturbations 
have three important advantages over the traditional 
Newtonian-gauge formalisms in terms of the proper mea- 
sures, e.g. [10, 17]. First of all, by Sec. II, the evolution of 
the coordinate, rather than proper, measures of pertur- 
bation is determined entirely by the physics within the 
local Hubble volume. 

The second, related to the first, improvement of the 
formalism is the explicit Cauchy structure of the new 
dynamical equations in the Newtonian gauge. Indeed, 
the traditional equations for the rate of change of proper 
density or intensity perturbations inevitably contain the 
time derivative of \& [10, 15-17, 39]. Yet, by the general- 
ized Poisson equation [c.f. eq. (40)], ^> itself depends on 
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the rate of change of the density and velocity perturba- 
tions of all the species. The terms or $ do not appear 
in the proposed equations. 

And the third benefit of quantifying density perturba- 
tion by d = 6n coo /n coo is a nonsingular generalized Pois- 
son equation for eq. (36). In terms of the traditional 
variables, 

V 2 * = 7 + 3 ^ ■ (4°) 

As t — > 0, the factor 7, eq. (34), diverges as H 2 ~ 1/t 2 . 
For the typical, e.g. inflationary, initial conditions this 
divergence on supcrhorizon scales is counterbalanced by 
a cancelation + 3Hu = 0(r 2 ). The seemingly singu- 
lar right-hand side of traditional eq. (40) may misguide 
an analytical analysis as well as cause significant numeri- 
cal errors and instabilities due to incomplete cancelation 
of the singular contributions. The regularity of in the 
supcrhorizon limit, on the other hand, is natural from 
eq. (36) for any conserved overdensities d a - There is gen- 
erally no large cancelation in the supcrhorizon limit of 
eq. (36). 

We now return to the question of whether the Newto- 
nian gauge provides a nonsingular description of super- 
horizon perturbations. As proved in Rcf. [40] , for any rea- 
sonable cosmology there exist perturbation solutions for 
which species' overdensities, the product Ti.u, and gravi- 
tational potentials in the Newtonian gauge remain finite 
in the supcrhorizon limit. In addition, eq. (31) admits 
singular solutions with u ~ 1/a. One might discard the 
singular modes on the grounds that if they are compa- 
rable to the regular ones when the perturbations exit 
the horizon then the singular modes become negligible 
soon after the exit. However, the solutions which ap- 
pear singular in the Newtonian gauge could be produced 
by a physical mechanism whose description is singular in 
the Newtonian but regular in a different gauge. Then 
on supcrhorizon scales in the Newtonian gauge decaying 
singular modes would appear to dominate frozen regular 
modes. 

This question cannot be resolved without taking into 
account dynamical equations. For example, we can imag- 
ine regular perturbations in metric whose shear does not 
vanish (V 2 e -/-> 0) in the supcrhorizon limit. The trans- 
formation of metric shear 

V 2 e~ = V 2 e + Sx l A (41) 

under a scalar gauge transformation shows that this met- 
ric requires a singular transformation to the Newtonian 
gauge, in which e ( Newt ) = 0. The corresponding New- 
tonian perturbations may then turn out singular. For- 
tunately, we know that, at least, slow-roll inflationary 
scenarios lead to superhorizon perturbations which are 
nonsingular in the Newtonian gauge [9]. Therefore, given 
the convenience of the Newtonian gauge, we use it for our 
analysis of post-inflationary cosmological evolution. 



IV. OVERVIEW OF PERTURBATION 
EVOLUTION 

A. Superhorizon evolution 

We can identify three popular descriptions of pertur- 
bation evolution on supcrhorizon scales. Following tra- 
dition, the general superhorizon perturbation is often 
viewed as a superposition of adiabatic and isocurvaturc 
modes, e.g. [41]. Adiabatic (also "curvature" or "isen- 
tropic") perturbations can be defined as perturbations 
for which the internal (proper) characteristics of all the 
matter species are uniform on some spatial slice. On this 
slice the metric remains perturbed. Isocurvature pertur- 
bations were idealistically considered to be the perturba- 
tions of the relative abundances of species while the total 
energy density was assumed unperturbed and the geom- 
etry homogeneous and isotropic. It has been, however, 
realized for two decades [42-46] that the decomposition 
into adiabatic and isocurvature modes generally is not 
invariant under cosmological evolution. Although an adi- 
abatic perturbation preserves its nature on superhorizon 
scales [4, 40, 47], a perturbation which appears isocurva- 
ture at a particular time, generally fails to be such in just 
a single e-folding of Hubble expansion. Despite this se- 
rious inconsistency, the picture of "isocurvature modes" 
continues to be influential, at least, in terminology. 

A different view of superhorizon evolution is given by 
the "separate universe" picture, e.g. [34, 48, 49]. In 
it, each Hubble volume is considered as an independent 
FRW universe, characterized by its own expansion rate 
and curvature. This description is self-consistent, and it 
explicitly exhibits the locality of superhorizon evolution. 
Nevertheless, it requires to keep track of "how far" each 
region has evolved in a global frame. The procedure for 
the latter is well defined but its outcome is not intuitive 
when the perturbations are not adiabatic. 

More recently, Wand et al. [20] noted that supcrhori- 
zon evolution of several uncoupled fluids becomes triv- 
ial if inhomogencity of a given fluid a is described by 
the perturbation of spatial curvature on the hypersur- 
faces of constant p a - Neglecting gravitational decays of 
species [21], these curvature perturbations, £ a , are indi- 
vidually conserved for uncoupled fluids even if the overall 
perturbation is not adiabatic. This approach lead to sig- 
nificant progress, e.g. [50, 51], in quantitative analysis of 
the curvaton-type scenarios [42-46]. 

The formalism presented in Sees. II and III further fa- 
cilitates the description of linear supcrhorizon evolution 
of multiple species as the common expansion of grav- 
itationally decoupled inhomogencitics in the individual 
species. It shows that if the inhomogeneities are quanti- 
fied by the perturbations of the coordinate number den- 
sities then their supcrhorizon freezing is explicit in any 
regular gauge. The conserved curvatures ( a [20] are sim- 
ply related to the coordinate density perturbations d a in 
the Newtonian gauge, as £ a = | e?i Ncwt ' ) [comment a) after 
eq. (29)]. The description in terms of d a does not require 
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to specify the conserved perturbations of each fluid a on 
a different set of hypersurfaces. This description retains 
its usefulness on all scales, including subhorizon. 

If the supcrhorizon cosmological perturbations are adi- 
abatic, as predicted by single-field inflation and favored 
by observations [52, 53], then the number density per- 
turbations d a of all the species are equal. This can be 
seen by evaluating the right-hand side of eq. (28) in the 
gauge in which all 5p a vanish. The common value of d a 
then equals three times the perturbation of the uniform- 
density slice curvature D, known as the Bardeen curva- 
ture C [6, 18]: d a = 3C- In other words, on superhorizon 
scales an arbitrary collection of adiabatically perturbed 
species is equivalent to a single fluid, having an equa- 
tion of state p(p) = p(p) and frozen coordinate number 
density perturbation d = 3£- 

Even if the overall cosmological perturbations are not 
adiabatic, any species which do not interact with the 
other species and are perturbed internally adiabatically 
(i.e., are unperturbed on a certain spatial slice, on which 
the other species may be perturbed) can be regarded as a 
single fluid, "same" in various regions beyond causal con- 
tact. Its density perturbation d a is also frozen (Sec. II C). 

The bulk peculiar velocities of species change with 
time even on supcrhorizon scales. In general, the evo- 
lution of velocities differs for various species. However, 
if supcrhorizon perturbations are adiabatic then the last 
two terms in eq. (31) vanish. (These terms are gauge- 
invariant and they apparently vanish for adiabatic per- 
turbations in the uniform density gauge.) Then for a 
Fourier perturbation mode with a wavevector k 

^-(au a )=a[<i> + 0(k 2 /H 2 )}. (42) 

Thus for adiabatic conditions the superhorizon accelera- 
tion of all the species is equal. 

B. Horizon entry 

When the perturbation scale becomes comparable to 
the Hubble horizon (the perturbations "enter the hori- 
zon"), the contribution of velocity divergence to d a , 
eq. (19), becomes appreciable. The evolution of the ve- 
locities is affected by metric inhomogeneities and the lat- 
ter are sourced by all the species which contribute to 
8T^ U . Thus, during horizon entry all such species influ- 
ence the evolution of all the density and velocity pertur- 
bations. 

After the entry, the perturbations of given species 
may or may not continue affecting other species gravi- 
tationally. The criterion for perturbations to continue 
being gravitationally relevant is evident from eqs. (36)- 
(37). On subhorizon scales these equations reduce to 
the Poisson equation fc 2 $ ~ k 2 ^ ~ — AnGa 2 (p + p)d = 
—4:irGa 2 J2 a (pa + p a )d a - If for some species after the 
horizon entry a 2 (p a +p a )d a decays, such as for relativis- 
tic species, then the species contribute negligibly to the 



gravitational potential. The gravitational impact of their 
perturbations can be ignored. On the other hand, if 
fl2 (Pa + Pa)d a ~ a 2 6 p a remains non- vanishing, such as 
for CDM during the matter era, then these species source 
metric perturbations and maintain their influence on the 
evolution of all cosmological species. 

V. FEATURES IN THE ANGULAR SPECTRUM 
OF THE CMB 

A. Preliminaries 

Prior to decoupling of photons from baryons at z ~ 
1100, Fourier modes of photon overdensity evolve accord- 
ing to wave equation (A12) in Appendix A. We find it 
useful to rewrite this equation as 

d\ + -^d\ + c 2 s k 2 (d~ / -B)=0, (43) 

where Rb = 3pb/(4p 7 ), c 2 = [3(1 +i?t,)] ~ 1 gives the sound 
speed in the photon-baryon plasma, 

B(r,k) = -3($ + * + i? b $), (44) 

and we continue to use the Newtonian gauge (26). Ac- 
cording to eq. (43), in each mode the density per- 
turbation d-y oscillates about its varying equilibrium 
value B(r,k). 

The presently observed CMB perturbations are given 
by a line-of-sight solution [15] of the radiation transport 
equation (33) complemented by the Thomson collision 
terms [10, 14, 54], eq. (A13). Denoting the collision terms 
by Ct, we have 

Ly + mVih = -4n,Vi($ + *) + C T . (45) 

Once the evolution of the dynamical variables cL, d v , d c , 
« 7 , etc. and the induced gravitational potentials $ and 
are determined, the line-of-sight integration of the radi- 
ation transport equation (45) becomes more straightfor- 
ward if the spatial derivatives of the potentials in eq. (45) 
are traded for their time derivatives as 

irff+niVitrff = 4($ + *)+C r Tj (46) 

where 

L cfi = i 7 + 4($ + *). (47) 

The variables i 7 (n) and t e ff(n) have identical multi- 
poles for all £ > 1. However, their monopolcs differ: 
(On = 1^7' ec l- (16), while (i e ff)n = 48 c ff, where the ef- 
fective temperature perturbation O c ff [55-57] in our vari- 
ables reads 

e cff = i d 7 + $ + (48) 

(In the traditional approach [56, 57], 6 off = ST^ Ncwt yT+ 
$.) Physically, O e fj quantifies the present temperature 
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perturbation of photon radiation which decoupled from 
a stationary thermal region with photon overdensity g? 7 
and potentials $ and 'S, provided the potentials subse- 
quently decayed adiabatically slowly over a time interval 
At A = k^ 1 . Indeed, under the described conditions 
i ff is conserved along the line of sight by eq. (46). 8 

The appearance of the term <f> + rather than the 
traditional <I>, in eq. (48) should not be surprising. The 
same sum $ + <F is also responsible for the integrated 
Sachs- Wolfe [58] (ISW) corrections to the intensity in 
the variable potential and for the deflection of photons 
by gravitational lensing. 9 

The full expression for the CMB temperature pertur- 
bation observed at the present time tq reads: 
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(49) 




kS(x)/jt with S(x)=x/3' 



FIG. 1: The evolution of a Fourier mode of CMB density per- 
turbation in the radiation era in terms of different measures: 
Coordinate number density perturbation d 7 = 8n~ t coo / n~ t coo 
(solid), proper density perturbation <5p 7 /p 7 in the Newtonian 
gauge (dashed), and effective temperature perturbation O e ff , 
eq. (48), (green solid), all normalized to coincide on small 
scales. 



B. Radiation era 



It follows from the line-of-sight integration [15] of eq. (46) 
with the scattering terms of eq. (A13). Here, g(r) is 
the probability for a photon emitted at time r to reach 
the observer unscattered, eq. (B5), and for the scalar 
perturbations Dj? = (VjVj — |<5yV 2 )7r p , where n p for 
Thomson scattering is given by eq. (A20). The pertur- 
bations ff, Vb, and ir p are evaluated along the line of 
sight t*(t) = (t — To)n, taking that the direction of AT 
observation is — n and at the observer's position r = 0. 
The complete formulas for the angular power spectra C; 
of CMB temperature and polarization induced by the 
scalar perturbations are summarized in Appendix B. 

We now proceed to connecting features in the angular 
spectrum of CMB temperature anisotropy to the internal 
nature of cosmological species. In Sec. VB we address 
the scales which enter the horizon during radiation dom- 
ination, and in Sec. V C the larger scales which enter 
during matter and dark energy domination. 



The variables (, e fj and © e ff are useful for constructing the CMB 
linc-of-sight transfer functions [15] from known solutions for per- 
turbation modes. Yet, these variables are not convenient as pri- 
mary dynamical variables and should be avoided when studying 
perturbations' gravitational interaction. The evolution equations 
for L e ff and © c g do not have an explicit Cauchy form and their 
dynamics does not satisfy condition I of Sec. I. 
Of the numerous existing formalisms for CMB dynamics in the 
Newtonian gauge, I am aware of only one [59, 60], by R. Durrcr, 
in that the impact of metric perturbations on radiation tempera- 
ture is described by <I> + I I', as opposed to <£>. Durrer's variable for 
linear perturbation of radiation intensity coincides with the lin- 
ear limit of our intensity perturbation t, constructed in Sec. II B 
for full theory. Yet, unlike the considered d a , her density per- 
turbations D a = (1 + w a )d a Newt ^ do not freeze for superhorizon 
perturbations. 



1. Radiation driving? 

It is often stated, e.g. [17, 25-29], that the CMB per- 
turbations on the scales of a degree and smaller (I > 200) 
arc significantly enhanced during the horizon entry in the 
radiation era. In a model without neutrinos, the claimed 
enhancement ( "radiation driving" ) of CMB temperature 
perturbations is by a factor of 3, translating into a fac- 
tor of 9 enhancement of the temperature autocorrelation 
spectrum. This conclusion is based on the apparent en- 
hancements of both o ff and ST/T in the Newtonian 
gauge during the entry (the solid green and dashed curves 
in Fig. 1). The enhancements are ascribed to a "resonant 
boost" of CMB perturbations by their own gravitational 
potential, whose decay soon after the horizon entry is 
said to be "timed to leave the photon fluid maximally 
compressed" [28]. 

However, in the concordance cosmological model as 
much as 41% of the energy density in the radiation era 
is carried by decoupled neutrinos. The evolution of free- 
streaming neutrinos and their contribution to metric in- 
homogeneities differs qualitatively [4] from the acoustic 
dynamics of the photon fluid. If without neutrinos the 
CMB perturbations were resonantly boosted by metric 
inhomogeneitics, neutrinos should be expected to sig- 
nificantly untune the "resonance" . In reality, neutri- 
nos decrease subhorizon ST/T by as little as 10%, as 
compared to subhorizon ST/T in a model with identical 
inflation-generated primordial perturbations but no neu- 
trinos. Moreover, the addition of noticeable density of 
early tracking quintessence [31, 61-63], whose perturba- 
tions evolve very dissimilar to the acoustic CMB pertur- 
bations (dashed vs solid lines in Fig. 2), not only fails to 
destroy the "resonance" but instead enhances subhori- 
zon 8T/T (Fig. 3 a). 
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We argue that the apparent radiation driving of photon 
temperature in the Newtonian gauge is a gauge artifact. 
First, the driving is absent when CMB density perturba- 
tions are quantified by the measure d 7 , satisfying condi- 
tions I and II. Indeed, in a photon-dominated universe 
the overdensity d 7 in the Newtonian gauge evolves as [4] 

/ 2sm<i _ CQS(p j ^ ip=kr/V3, (50) 

as follows from eqs. (All) and (36). Hence, without 
neutrinos or early quintessence, the amplitude of sub- 
horizon oscillations of photon overdensity, Amp[d 7 ], ex- 
actly equals the primordial photon overdensity, d 7 (r — ► 
0) = c?m- By Sec. II D, any other measure of photon per- 
turbations d 7 that satisfies conditions I and II coincides 
with dry on super- and subhorizon scales. Any such mea- 
sure, therefore, would also detect no radiation driving: 
Amp[d 7 ] =d;( T ->0). 10 

A formalism-independent distinction between physical 
driving of perturbations and a gauge artifact can be made 
by comparing the studied model to a reference model in 
which the metric becomes homogeneous by reentry. In 
the latter model any driving by metric inhomogencities 
at reentry is objectively absent. A consistent procedure 
to achieve homogeneous geometry before reentry for any 
given primordial perturbations of arbitrary species was 
described in Sec. II D. In the homogeneous (FRW) metric 
the evolution of photon fluid overdensity readily follows 
from eq. (A12) with c\ = 1/3 as 



j(undriven) j 

dry — tijr 



COS (f. 



(51) 



Thus if in the superhorizon limit the photon and metric 
perturbations in the photon-dominated and "undriven" 
models are identical, after horizon entry the acoustic os- 
cillations in the two models have equal amplitude. There 
is no causal mechanism for a gravitational boost of the 
amplitude in the second model. Therefore, the ampli- 
tude should also not be considered boosted in the photon- 
dominated model. 

The 3-fold increase of the Newtonian-gauge perturba- 
tion STj ewt yT 7 = | dry + ^ in the second of the above 
models occurs when the metric becomes homogeneous 
on superhorizon scales. This increase can hardly be at- 
tributed to a resonant boost of the acoustic oscillations 
from the decay of their gravitational potential at the 
maximum photon compression, as stated in [28]. The 

enhancement of the perturbation measure is 
possible because its evolution is not determined only by 
local dynamics. 



10 For comparison, the perturbation of proper density in the New- 
tonian gauge is 



6Ty Ncwt ^ ( 2 sin ip 2 cos ip 2 sin <p 

■ Tf, = d in COS ip -f- — 

Try \ <fi ^ l lp J 




k S(t) / n 



FIG. 2: a) Density perturbations d a of typical cosmological 
species in the radiation epoch dominated by tightly coupled 
photons (7, red solid) and 3 flavors of decoupled relativistic 
neutrinos (y, dash-dotted). The subdominant species are cold 
dark matter (CDM, black solid), and a scalar field (tracking 
quintessence Q, dashed). The primordial perturbations are 
assumed adiabatic. The dotted curve shows the equilibrium 
value of d 1 oscillations, eq. (44). b) Bulk peculiar velocities 
of the same species. 



Gravitational impact 0} perturbations on the CMB and 
CDM 



Fig. 2 a) displays the linear evolution of coordinate 
density perturbations d a (r, k) in a photon-baryon fluid, 
decoupled relativistic neutrinos, CDM, and a mini- 
mally coupled scalar field [tracking quintessence, Q, with 
w q( z ) = ^totai(z)]- Fig. 2b) shows the corresponding 
bulk velocities of the species. The plots describe the ra- 
diation epoch, dominated by photons and three standard 
flavors of neutrino. They are obtained by the numer- 
ical integration of eqs. (A1)-(A26) and (36)-(37) with 
the regular adiabatic initial conditions, normalized to 
d(k/H — ► 0) = 1. The chosen evolution parameter is a 
dimcnsionlcss phase of photon acoustic oscillations kS(r) 
where S(t) = C c s dr = r/y/3 (for Rb <C 1) is the size of 
the acoustic horizon. 

In the description of the acoustic dynamics of the 
CMB by eq. (43), the gravitational impact of the species' 
perturbations is mediated through the "restoring force" 
-c 2 s k 2 (dry - B), with 



B(r,k) = -3($ + * + 



(52) 



For it, Arnp[<5n 7 /n 7 ] = 3 lim T ^o <5n 7 /n 7 



The potentials $ and \& are sourced by perturbations of 
densities, peculiar velocities, and anisotropic stress of the 
species as described by eqs. (36)- (37). 

For typical cosmological species, the growth of sub- 
horizon perturbations during the radiation era is insuffi- 
cient for supporting a non-vanishing gravitational poten- 
tial. Hence, the equilibrium photon overdensity B(t, k) 
on subhorizon scales tends to zero. Thus in the radiation 
era the gravitational impact of perturbations is confined 
to horizon entry. As for the later subhorizon CMB oscil- 
lations, this impact can only rescale the oscillations' am- 
plitude and additively shift the phase. We now consider 
the gravitational impact on the CMB that is caused by 
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decoupled rclativistic neutrinos and a scalar field (early 
quintessence) . 

Addition of extra ultrarelativistic neutrinos changes B 
as shown by the solid curve in Fig. 3 c). The increase 
of N v reduces B during the initial growth of d 1 while 
enhances B for a significant period of d 1 fall. The re- 
sulting modulation of the restoring force —c\k 2 [d 1 — B) 
somewhat damps d 1 oscillation. The reduction of the 
amplitude of CMB oscillations by neutrino perturba- 
tions [4, 25] is evident from the thin solid curve in 
Fig. 3 a). 

On the contrary, the density perturbation of canonical 
tracking quintessence [31, 61-63] is significantly reduced 
right after horizon entry [see the dashed line in Fig. 3 a) .] 
This reduction causes a prominent drop in B(r), shown 
by the dashed curve in Fig. 3c), at the time when d a is 
on average negative. As a result, the acoustic oscillations 
gain energy and their amplitude increases [64]. 

The qualitative difference in the evolution of density 
perturbations for quintessence versus the other consid- 
ered species can be understood as follows. For the regular 
adiabatic initial conditions, the number density pertur- 
bations d a of all the species are equal and the velocities 
vanish beyond the horizon (kS -C 1). Moreover, the 
initial accelerations of all the species, eq. (42), are also 
equal, as confirmed by Fig. 2 b). However, the rate of 
change of d a at reentry is very sensitive to whether or 
not pressure of the species is fully specified by their local 
density. If it is, such as for photons, CDM, and neu- 
trinos, the "nonadiabatic" term on the right-hand side 
of eq. (19) vanishes. The corresponding <i a 's agree in 
their evolution up to 0(/c 4 t 4 ) corrections. On the other 
hand, for a classical field the nonadiabatic term becomes 
non-zero during the entry [c.f. the first of eqs. (A26).] 
As a consequence, quintessence density perturbation dQ 
[dashed line in Fig. 2 a)] noticeably deviates from the per- 
turbations of the other species already in 0(fc 2 r 2 ) order. 

Quantitatively, a change of N v from 3 to 4 changes 
the amplitude of d 1 oscillations by —1.7%. Addition of 
equal energy in the form of tracking quintessence, on the 
contrary, enhances the amplitude by 1.6%. 

Perturbations of CDM evolve according to eqs. (Al) 

as 

d c = - J (Vi<) dr, avl = - J (aVj$) dr. (53) 

The impact of metric perturbations on CDM is thus de- 
termined by a normalized quantity — (a/aent)^, where 
a cnt = a\k T =i- Fig. 3d) shows how this quantity 
is affected by the addition of extra neutrinos and 
quintessence. By comparing Fig. 3 b) and Fig. 3 a), we 
see that the signs of the gravitational impact of pertur- 
bations in either neutrinos [4, 25] or quintessence [65] on 
the power of CDM are opposite to the signs of the impact 
of these species on the acoustic amplitude of the CMB. 
Thus on the scales entering the horizon before matter- 
radiation equality (£ > 200, k > 10 _2 Mpc~ 1 ) the ratio 



of CMB to CDM powers, which is insensitive to the un- 
known primordial power of cosmological perturbations, is 
slightly decreased by decoupled particles. On the other 
hand, this ratio is increased by tracking quintessence. 

The gravitational impact of perturbations on the phase 
of CMB oscillations is easier to analyze in real space [66- 
69]. The phase shift, along with all the other character- 
istics of subhorizon oscillations, is fully specified by the 
singularity of the CMB Green's function at the acous- 
tic horizon [4]. An analysis of this singularity [4] shows 
that for adiabatic initial conditions the phase is affected 
only by the perturbations which propagate faster than 
the acoustic speed of the photon fluid. 

Perturbations of both ultrarelativistic neutrinos and 
quintessence propagate at the speed of light c, i.e., faster 
than the acoustic speed c s < c/V3- The phase shift 
induced by neutrinos displaces all the CMB peaks in 
the temperature and E-polarization spectra by an about 
equal multipole number [4] 

Al„ « -57 Pv . (54) 

Pj + Pv 

The three standard types of neutrino give « —23, 
large enough to be observable with the present data [70- 
76]. One additional fully populated type of relativistic 
neutrino increases the shift by [4] A/+„ ~ —3.4. Accord- 
ing to numerical integration of eqs. (A1)-(A26) and (36)- 
(37), the same additional radiation density (AA c ff = +1) 
in tracking quintessence leads to an almost triple shift 
AIq ~ -11. 

While in a formal treatment of Boltzmann hierar- 
chy for the perturbations of free-streaming neutrinos, 
e.g. [77], the anisotropic stress might appear to be the pri- 
mary cause of their characteristic gravitational imprints, 
this conclusion is misleading. For example, anisotropic 
stress of minimally coupled tracking quintessence, whose 
perturbations also change the amplitude and phase of 
acoustic oscillations, is identically zero. In particular, 
a real space approach [4] shows that for adiabatic per- 
turbations the phase is shifted if and only if some per- 
turbations physically propagate faster than the acoustic 
speed c s . 

C. Matter era 

1. 25-fold Sachs-Wolfe suppression 

In the epoch dominated by pressureless CDM and 
baryons which have decoupled from the CMB, metric 
perturbations (gravitational potentials) do not decay on 
subhorizon scales. Through the finite potentials, matter 
continues to influence CMB long after horizon entry. For 
the perturbations entering the horizon in this epoch, mat- 
ter overdensity grows linearly in the scale factor a oc t 2 , 
as d = dj n (l + £; 2 t 2 /30). If the initial conditions are 
adiabatic, photon and neutrino perturbations contribute 
negligibly to inhomogeneitics of the metric in the matter 
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FIG. 3: RADIATION ERA. The impact of perturbations in 
decoupled neutrinos and quintessence on the CMB and CDM. 
The wide solid curves give d 7 (a) and dcoM (b) in the epoch 
dominated by tightly coupled photons and three standard fla- 
vors of decoupled neutrinos. The thin solid curves describe 
a model with doubled neutrino energy density. The dashed 
curves describe a model with the standard neutrino content 
plus the same energy density in tracking quintessence, with 
wq = 1/3. The bottom panels demonstrate how these modi- 
fications of the fiducial model change the gravitational terms 
which drive the CMB (c) and CDM (d). 



era. Then Newtonian gravitational potentials, induced 
entirely by perturbations of nonrelativistic matter, re- 
main time-independent throughout all the linear evolu- 
tion: 



* = * = ~dm, 
5 



(55) 



as can be verified with eqs. (36) and (Al). 

The acoustic dynamics and subsequent free-streaming 
of CMB photons in the gravitational wells of non- 
relativistic matter are well understood on subhorizon 
scales, e.g. [17, 25, 26, 28, 29, 78]. However, the tra- 
ditional description of horizon entry in the matter era, 
e.g. [26, 28, 39, 57], requires corrections. 

For the considered perturbations, entering the horizon 
after CMB last scattering but long before dark energy 
domination, the Doppler and polarization contributions 
to the observed CMB anisotropy AT/T [second and third 
terms in eq. (49)] are negligible. The ISW contribution 
>i>+ v I' is also absent. Indeed, the potentials are static dur- 
ing the matter era and their later decay due to the dark 
energy domination is for these modes adiabatically slow. 
The remaining contribution to the observed AT/T (49) 
is then described by the effective temperature perturba- 
tion e cfr (48), 



AT ST 
— - Bcff - -jr 



2$. 



(56) 



FIG. 4: MATTER ERA. The perturbations of the CMB den- 
sity d 1 and effective temperature O c fi (a) and of CDM den- 
sity d c (b) in the matter era dominated by pressureless mat- 
ter (solid curves) and by equal densities of pressureless mat- 
ter and tracking quintessence, with wq = (dashed curves). 
Similarly to Fig. 3, the bottom panels show the change in- 
duced by quintessence in the gravitational terms which drive 
the CMB (c) and CDM (d). 



Here, wc defined 



ST 
~T 



— g d-y , in ; 



(57) 



with <2 7) i n = dj(k <C H(r)) and set $ = in the matter 
era. The quantity (57) is a natural measure of the in- 
trinsic CMB temperature perturbation on superhorizon 
scales. This is the perturbation to be always observed 
for the decoupled CMB photons if the metric becomes 
homogeneous on superhorizon scales and remains such 
during the subsequent evolution. 

If the initial conditions are adiabatic (Sachs- Wolfe 
case) then, by eq. (55), $ = — gd in = — gd^^. Thus 
the classical Sachs- Wolfe result AT/T = $/3 [58] is due 
to a partial cancellation of the primordial CMB density 
perturbation ST/T\- m = — 1$ against the gravitational 
redshift 2$ in proportion —5 : 6. The 5-fold suppression 
of the primordial temperature perturbation results in an 
astonishing (and experimentally evident) suppression of 
the angular power of CMB temperature Ci by a factor 
of 25. 

The traditional interpretations of the Sachs- Wolfe re- 
sult also note certain, 2-fold rather than 5-fold, suppres- 
sion of the intrinsic CMB temperature fluctuations by the 
matter potential. The relation AT/T = $/3 is usually 
interpreted as a partial cancellation between the "intrin- 



sic" temperature anisotropy <5T( Nowt )/T|i n 
the gravitational redshift $ [26, 28, 39, 57]: 



r$ and 



AT 
~T~ 



Jji(Newt) 



$ (misleading!) (58) 
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Let us present two examples where such quasi- Newtonian 
considerations lead to wrong expectations. 

One example is a classical scenario of "isocurvature" 
initial conditions, e.g., from a phase transition early in 
the radiation era, when p 7 3> p c . If immediately after the 
phase transition 5p 1 + 5p c = then |5p 7 /p 7 | <C \Sp c /p c \. 
For a fixed magnitude of cosmological inhomogencities 
\$Pc/pc\ ~ 10 -5 , the primordial photon perturbations 
Sp-y/pj can then be ignored. However, it would be 
wrong to conclude that (5T( Ncwt )/T| in = and that from 
eq. (58), due to the gravitational redshift in the matter 
era, AT/T| isoc = $ (incorrect!). The correct answer is 
AT/T| isoc = 2$, in agreement with eq. (56) where now 
6T/T\. m = 0. 

Another example where the interpretation (58) mis- 
leads is the opposite case of non-zero 5p 7 /p 7 but homo- 
geneous and isotropic geometry in the matter era. If 
the primordial perturbations of radiation in this scenario 
are identical to the perturbations in the adiabatic model 
then the interpretation (58) suggests that now, when 
the "2-fold" suppression by 4> is absent, AT/T| $=0 = 
—2 AT/T| adiab (incorrect!). The correct answer, imme- 
diate from eqs. (56) and (57), is 



AT 
~~T~ 



ST 
~~T 



AT 



(59) 



adiab 



in agreement with our claim of 5-fold suppression of 
AT/T by $ in the adiabatic case. 

The homogeneous geometry of the last scenario could 
be generated with the matter initial conditions d m , in ~ 0, 
yet there is no physical reason to expect such initial con- 
ditions [38]. However, d c and the matter-induced po- 
tentials $ and \& may be partly smoothed early during 
the horizon entry by non-standard dynamics of the dark 
sector. In the limiting case of a perfectly homogeneous 
metric since very early time, we recover the unsuppressed 
CMB oscillations, with the amplitude d 7i i n in overdensity 
and AT/T = ST/T\ in = -5 AT/T| adiab in temperature. 



2. Implications 

The observed CMB multipoles with I < 200 are pri- 
marily generated by the perturbation modes which enter 
the horizon during the matter era. These modes evolve 
in a static potential of collapsing CDM and baryons. The 
potential suppresses the corresponding C/'s by an order 
of magnitude, see the solid curve in Fig. 5 for the "con- 
cordance" ACDM model [76, 79-81]. 

The suppression is somewhat reduced for I < 10 be- 
cause of the modes which enter the horizon when dark 
energy becomes dominant and the potential decays. For 
these modes the enhancement of CMB power is fully de- 
scribed by the ISW term $ + ^, eq. (49). In general, 
however, the reduction of the 25-fold Sachs- Wolfe sup- 
pression is not equivalent to the ISW enhancement. An 
extreme example is the scenario in which $ + ^ vanishes 
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FIG. 5: Suppression of the angular power spectrum of CMB 
temperature Ci for I < 200 by dark matter perturbations. 
The figure compares three models which differ only by per- 
turbations in the dark sector: The "concordance" ACDM 
model [76, 79-81] (solid); a model where CDM and cosmolog- 
ical constant are replaced by a dynamical scalar field with the 
same background p(z) (dashed, the overall power of pertur- 
bations is decreased to match the concordance C300); finally, 
a ACDM model where CDM density is artificially arranged 
to be unperturbed beyond the horizon (dash-dotted, the pri- 
mordial power of photon, baryon, and neutrino perturbations 
is the same as in the concordance model). 



prior to CMB decoupling. Given equal primordial per- 
turbations, the observed CMB power C; in this scenario 
exceeds the Sachs- Wolfe power by 25 times, although in 
both scenarios the ISW effect is absent. 

More realistically, the density perturbation d a of some 
of the dominant species in the matter era may grow 
slower than [^(p+p)}- 1 ~ 1/H 2 ~ t 2 . Then, by eq. (36), 
these species contribute only decaying terms to the grav- 
itational potentials. This diminishes the suppression of 
the CMB power at I < 200. As an illustration, the upper 
solid and dashed curves in Fig. 4 a) show d 1 in a CDM- 
dominated universe (Sachs- Wolfe case) and in a model 
where the energy density of the dark sector is split equally 
between CDM and quintessence whose background den- 
sity tracks CDM. In the second model the CMB tem- 
perature perturbations on subhorizon scales are larger 
by 93%. 11 

Thus any physics which reduces metric inhomo- 
geneities in the matter era strongly influences CMB tem- 
perature anisotropy. Some of the influence may be due 
to a different ISW effect from the temporal change of the 
gravitational potentials. It can also be due to reduced 
gravitational lensing of the CMB, appearing in the sec- 
ond perturbative order. However, the most prominent 
impact is caused by the very reduction of the potentials, 
suppressing the auto-correlation of CMB temperature by 
a factor of 25. 



11 If in the matter era some of the dominant species do not cluster, 
the gravity-driven growth of CDM inhomogeneities also signifi- 
cantly weakens. This is seen distinctly in Fig. 4 b). A classical 
example of such species is baryons prior to their decoupling from 
the CMB, e.g. Ref. [25]. 
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Two examples are shown in Fig. 5. The figure displays 
the CMB temperature power spectra C\ obtained with a 
modified CMBFAST code [15] for three models which differ 
only by the dynamics of perturbations in the dark sec- 
tors. The solid curve describes the concordance ACDM 
model [76, 79-81]. The dashed curve refers to a model 
where the dark sector is composed entirely of a classical 
canonical scalar field Q. In this model at every rcdshift 
the field background density Pq(z) is chosen to equal the 
combined CDM and vacuum densities p c (z) + pa of the 
ACDM model. However, since the field perturbations 
do not grow, the metric inhomogeneities of the second 
model decay after the horizon entry in the matter era. 
As a result, its CMB power for I < 200 is enhanced. 

The last, dash-dotted, curve in Fig. 5 describes a model 
which has the matter content of the concordance ACDM 
model. However, while in the previous two scenarios the 
initial conditions were adiabatic, here CDM density per- 
turbation d c is set to zero on superhorizon scales (the 
initial values for d 7 , d v , and db are unchanged). Then, 
again, the metric inhomogeneities in the matter era are 
reduced. While this model has the least ISW effect, the 
smoother metric causes the largest enhancement of the 
CMB power among the compared models for £ < 200. 



VI. SUMMARY 

There is a variety of descriptions of inhomogeneous 
cosmological evolution which are apparently dissimilar 
on the scales exceeding the Hubble scale. To single out 
the descriptions least afflicted by gauge artifacts, we sug- 
gest considering the measures of perturbations that meet 
the following criteria: I. The measure's dynamics is de- 
termined entirely by the physics within the local Hub- 
ble volume. II. The measure is practically applicable in 
Minkowski limit and retains its physical meaning in any 
FRW geometry on all scales. 

We showed (Sec. II A) that in the fully nonlinear the- 
ory these criteria are fulfilled for a perturbation of one- 
particle phase-space distribution /(x^, Pi), considered as 
a function of particle canonical momenta Pi. The of- 
ten used perturbations of f{x ll ,p) or f(x^ 7 q), where p 
is the proper momentum and q = ap, do not satisfy re- 
quirement I. Perturbed nonlinear dynamics of radiation 
can be appropriately described by the intensity pertur- 
bation t(x A ',ni), constructed by integration of f(x^,Pi), 
eqs. (7) and (9). Finally, a description of perturbed lin- 
ear dynamics of arbitrary species satisfies criteria I and II 
if the perturbations are quantified by coordinate number 
density perturbations d = Sn coo /n coo , eq. (20). In linear 
theory these measures of perturbations in phase-space 
distribution, intensity, and density are interrelated with 
eqs. (24) and (25). 

Similarly to the traditional proper measures of per- 
turbations, e.g. AT(h)/T and Sp/p, the proposed "co- 
ordinate" measures t(ni) and d = Sn coo /n coo are gauge- 
dependent. Yet the gauge dependence of the proposed, 



but not traditional, measures disappears in the super- 
horizon limit. More generally, we prove that on super- 
horizon scales the values of the density perturbations are 
the same in terms of any measure satisfying the formal 
versions of criteria I and II, namely, conditions I' and II' 
on page 7. 

In the Newtonian gauge (26) the suggested variables 
provide additional technical advantages. One is the 
explicit Cauchy structure of the dynamical equations 
(cq. (19), Sec. IIIB, and Appendix A). Another is a 
nonsingular superhorizon limit of the generalized grav- 
itational Poisson equation (Sec. IIIC). Analytical so- 
lutions in the new formulation inevitably appear more 
tractable [4]. As an example, compare the solutions for 
density perturbations in the photon-dominated model, 
eq. (50) versus the traditional solution in Footnote 10. 

In the proposed class of formalisms, classical linear 
density perturbations of individual species manifestly 
decouple gravitationally on superhorizon scales. This 
view of linear superhorizon evolution of multiple species 
(Sec. IV A) agrees with the picture provided by the con- 
served curvatures of Wands et al. [20] . The density per- 
turbations of the species recouple during the horizon en- 
try (Sec. IV B). 

On the other hand, our description reveals a different 
from the traditionally assumed origin of several promi- 
nent features in the angular power of CMB tempera- 
ture Ce- First of all, the apparent "resonant boost" of 
CMB perturbations for £ > 200, presumably by their 
self-gravity during the horizon entry in the radiation 
era [17, 25-29], is found to be a gauge artifact of the 
Newtonian gauge. We demonstrated (Sec.VB 1) that the 
driving is absent when density perturbations are quanti- 
fied by any measure satisfying conditions I' and II'. 

We showed (Sec. VB1) that the "radiation driving" 
picture leads to incorrect expectations for the observed 
CMB signatures of dark species. Since the resonant self- 
boost of CMB perturbations does not physically exist 
and cannot be untuned by the gravitational influence of 
non-photon species present in the radiation era with non- 
negligible densities, the CMB signatures of such species 
are not anomalously amplified. This corrected expecta- 
tion is confirmed by calculations for decoupled relativis- 
tic neutrinos and for tracking quintessence (Rcf. [4] and 
Sec. VB2). 

A moderate impact of such species on the CMB power 
at small scales, less affected by cosmic variance, neverthe- 
less already allows meaningful constraints on the species' 
abundance and nature, Refs. [70-77]. The constraints 
will significantly improve in the near future [4, 82, 83], 
when the smaller angular scales become resolved. 

Notably, the CMB is sensitive to the propagation ve- 
locity of small-scale inhomogeneities in the dark sec- 
tors [2]. For neutrinos, quintessence, and any other 
species for which this velocity exceeds the sound speed 
in the photon fluid c s w c/%/3, the perturbations cause a 
nondcgcncratc uniform shift of the acoustic peaks [4, 84] . 
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The CMB power is also sensitive to nonadiabatic pres- 
sure, generally developing for a classical field. As dis- 
cussed in Sec. VB2, under adiabatic initial conditions, 
the nonadiabatic pressure affects the evolution of species' 
overdensities already at the order 0(k 2 r 2 ) while the reg- 
ular pressure enters only at 0(fc 4 T 4 ). In the example 
of tracking quintessence, nonadiabatic pressure causes 
early suppression of quintessence density perturbations 
(Sec. VB2 and Fig. 2). This enhances the CMB power 
while neutrino perturbations reduce it (Fig. 3). 

Although CMB perturbations are not resonantly 
boosted on small scales, the CMB multipoles with £ < 
200 are significantly suppressed by the gravitational po- 
tential which does not decay after the horizon entry in the 
matter era. The contribution of the perturbations which 
enter the horizon when \Sp/Sp\ <C 1 and the potential is 
time- independent, to the angular power of CMB temper- 
ature is suppressed by 25 times (Sec. VC1). The sup- 
pression is thus much more prominent than a factor of 4 
reduction of the CMB power in the conventional inter- 
pretations [26, 28, 39, 57] of the Sachs- Wolfe effect [58]. 12 

As a consequence, any mechanism which smoothes 
metric inhomogeneity after radiation-matter equality at 
redshift z ~ 3500 significantly enhances CMB anisotropy 
(Sec. VC2). Examples of such mechanisms are a notice- 
able fraction of quintessence at all post-equality rcdshifts, 
interaction or unification of dark matter and dark en- 
ergy, e.g. [30, 86-89], or alternative explanations for the 
apparent signatures of nonbaryonic matter, e.g. MOND- 
inspired theory of gravity [90]. 13 The enhancement is 
accounted for by the ISW effect when the metric is 
smoothed after CMB last scattering at z ~ 1100. How- 
ever, even if the ISW effect is absent, but the metric at or 
before CMB last scattering is already smoother than in 
the standard scenario, the enhancement exists. Thus the 
experimentally evident suppression of CMB multipoles 
with £ < 200 severely constrains the mechanisms which 



12 The CDM-induced suppression of the CMB power at low £'s is 
complimentary to another characteristic impact of dark matter 
potential on the CMB known as "baryon drag" or "baryon load- 
ing" [25, 39]. The baryon loading is the enhancement of the 
odd acoustic peaks in Ci relative to the even ones due to the 
shift of the zero of © e ff oscillations by —R^, eqs. (48) and (44). 
The alternation of the peak heights is generally absent in CDM- 
free models, e.g. [85]. While the suppression of C; at low Ps 
is controled entirely by the metric inhomogeneity <E>, the baryon 
loading is, in addition, proportional to baryon density p b oc Q. b h 2 . 

13 After the first appearance of this work (astro-ph/0405157) the 
predicted enhancement of CMB power spectrum at low I's for 
the Bekenstein TeVeS [90] implementation of MOND [91] has 
been observed in a numerical analysis of Skordis et al. [92], The 
authors of this analysis, nevertheless, ascribe the enhancement 
entirely to the ISW effect. This interpretation is generally in- 
correct. If by CMB last scattering the metric is perturbed less 
than in the standard ACDM model, the ISW effect in the models 
with enhanced CMB power may be smaller than in the standard 
model. 



reduce metric perturbations at any time in the matter 
era. 
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APPENDIX A: PERTURB ATIVE DYNAMICS OF 
TYPICAL SPECIES 

In this Appendix we give the closed equations for the 
linear evolution of scalar perturbations of typical cos- 
mological species: cold dark matter (c), decoupled neu- 
trinos (v), photon-baryon plasma (76), and a classical 
scalar field, or quintessence, (<p). We use the Newtonian 
gauge (26), in which the gravitational potentials <E> and iff 
arc determined by the perturbations of the total density, 
velocity, and anisotropic stress with eqs. (36)-(37). 



1. Cold dark matter (CDM) 

Before CDM streams cross, in particular, in the lin- 
ear regime, CDM is equivalent to a fluid with vanishing 
pressure and vanishing anisotropic stress (p c = Sp c = 0, 
a c = 0). Then from eqs. (19) and (31) 

d c = V 2 u c , u c = ~Hu c + (Al) 

Elimination of u c gives a single second-order equation for 
CDM density perturbation 

d c + Ud c = V 2 $. (A2) 



2. Neutrinos 

a. General case 

When the universe cools below 1 MeV, weak interac- 
tions of the standard neutrinos become cosmologically ir- 
relevant. Then the full evolution of neutrino phase-space 
distribution / is governed by collisionless Boltzmann 
equation (2). Straightforward linearization of eq. (2) 
shows that the scalar component of / perturbation 

dfix^Pd^fix^Pil-fiP) (A3) 
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evolves in the linear order in the Newtonian gauge (26) 
as 



(df) + V t V, (df) = E J| V t V, ($ ■ 



(A4) 



where E = (P 2 + mV) 1 / 2 and V l = Pi/E. 

The initial conditions for the perturbation df on su- 
pcrhorizon scales can be determined as follows. Consider 
spatial hypersurfaces of constant radiation density p r 
and, to be specific, assume that neutrino chemical poten- 
tial is unperturbed. Then prior to neutrino decoupling 
the neutrino phase-space distribution / pr (p P r) in proper 
variables d 3 x pr d 3 p pr is spatially uniform (independent 
of x l ) on a hypersurface p r = const. In a gauge of uniform 
radiation density and zero shear, in which p r = const de- 
fines the hypersurfaces of constant time and the 3-metric 
is gij = a (l + 2f r )5y, the proper and canonical neutrino 
momenta are related as p pl = Pi/[a(\ + ( r )]. Thus in the 
variables (x z , Pi) 



f(x\Pi)=f pi 



P 



a[l + Cr(^)] 



(A5) 



(We used that d 3 x pr d 3 p pv = d 3 x l d 3 Pi.) The correspond- 
ing linear perturbation of fix 1 , Pi) is 



Note that substitution of the last result in cq. (24) or 
integration of eq. (A5) over d 3 Pi returns d v = 3£ r . 

Since f(x l 7 Pi) is a 4-scalar and for the decoupled par- 
ticles f(P) is time- independent, df is a gauge invariant. 
Therefore, eq. (A6) gives the initial conditions for df evo- 
lution in any gauge, including the Newtonian. 



Ultrarelativistic limit 



Perturbation of the intensity of ultrarelativistic neu- 
trinos is conveniently described by the variable t(x M ,7ii), 
introduced by eq. (9). After neutrino decoupling, the full 
and the linear evolution of neutrino intensity perturba- 
tion is given by eqs. (10) and (15) respectively. The lat- 
ter shows that for linear scalar modes in the Newtonian 
gauge 



Lv + riiVii v = -4n,V,($ + *). 



(A7) 



Scalar perturbations of i v (n,i) may be parameterized 
by scalar multipole potentials d„i as [4] 



\i{m) = D{m) 



2=0 



where Pi(z) are the Legendre polynomials. The po- 
tentials d v i with I — 0,1, and 2 describe the neu- 
trino number density perturbation, velocity potential, 



and anisotropic stress potential respectively: d v § = d v , 
d v \ = u v , and g?„ 2 = °~v [4]. In the special case of a 
plane-wave perturbation l(x) oc exp(ifc-:r), W l di are 
proportional to the scalar (m = 0) multipoles liq = 
J d 2 nY l Q(rii)L(ni), where the polar direction is taken 
along k. By eq. (A7) and the identity (21 + l)pPi(p) = 
(I + l)Pi+i(n) + IPi-i(p), the neutrino multipole poten- 
tials evolve as 

d, A = d v ,i-! + V 2 d„ ii+ i + + *) (A9) 

(Sn in the last term is the Kronecker symbol.) 



3. Photon-baryon plasma 

a. Tight coupling limit 

Perturbations in the tightly coupled photon-baryon 
plasma (u~ ( = Ub, er 7 & = 0) propagate as sound waves 
with an acoustic speed, e.g. [17], 



8p 7 b 
5p 



'yb 



p lb 3(1 + Rb)' 



(A10) 



Here Rb = Pb/(p-y + P 7 ) = 3pb/(4p 7 ) and the photon- 
baryon perturbations are assumed internally adiabatic 
(A6) (^7 — db)- In this limit, from eq. (19) and eqs. (31), (27), 



<i 7 = V 2 u 7 , 



c 2 d 7 - H(l - 3c 2 )w 7 + $ + 3c 2 *. 



(All) 



Elimination of u 7 yields 



d-y + H(l - 3c 2 ) dy - c 2 V 2 d 7 = V 2 ($ + 3c 2 *). (A12) 

To describe photon and baryon inhomogencities be- 
yond the tight coupling limit we should confront direc- 
tional variation of photon intensity and account for the 
dependence of photon-baryon scattering on photon po- 
larization. 



b. Photon intensity 

The scalar component of the perturbation t 7 (ni) of 
polarization-summed photon intensity evolves according 
to the transport equation (33) with additional Thomson 
scattering terms [10, 14, 54]: 



n,Vii 7 



1 



-4niVi($ + *) 
4 



(A13) 



3^ 



£(-!)' (2Z + 1)^(71^ ) V 'd ; , (A8) Here, 



(an e cr T1 - 



lomson 



{n^if 



(A14) 



is a mean conformal time of a photon collisionlcss flight, 
v' b = — ViUb is the velocity of baryons, and ir p is defined 
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below by cq. (A20). Similarly to the case of neutrinos, 
we introduce scalar multipoles d~i of photon intensity 
perturbation t 7 (n,) according to expansion (A8). Then 
transport equation (A13) expands into a hierarchy of 
multipole evolution 14 

ii-y = i d 1 + |V 2 <T 7 + $ + * - J- (uj - u b ) , 

2 3 V72 J 1 /" 1 \ (•^■1^) 

cr 7 = g w 7 + 5 V d 7 , 3 - — - io M » 

where for the last equation d 7 2 = cr 7 . 

c. Photon polarization 

The phase-space density of partially polarized pho- 
tons with spacetime position and momentum p is 
described by a 2 x 2 Hermitian tensor fapix^, Pi), whose 
indexes a and /3 refer to two spatial directions orthogonal 
to Pj. Given a polarimeter selecting the photons whose 
electric field varies as E cx ecxpiiPiX 1 ), where e is a 
complex unit 3- vector orthogonal to Pj (gije*' l e J = 1 and 
P^e* = 0), the tensor / a| g specifics the phase-space den- 
sity of the photons passing the polarimeter as e* a f a p£ ■ 

We define the intensity tensor of partially polarized 
photons by integrating fapix^ , Pi) over P 3 dP, similarly 
to eq. (7): 

I a0 (x", ni ) = / P^dPfc^^iP). (A16) 
Jo 

The connection of I a p to the energy-momentum tensor of 
photons with a specific polarization e is given by eq. (6) 
with / — * e* a f a 0e^. Finally, we generalize the intensity 
perturbation i of eq. (9) to polarized radiation as 

La p( x ^ ni ) = ^-1, (A17) 

where I = Tr I a p is the polarization-summed background 
intensity. 

We will restrict the discussion of polarization evolu- 
tion to harmonic plane- wave perturbations, which vary 
in space as exp (ik ■ x). While the general analysis is also 
possible, we defer it to a separate publication. In linear 
theory, solving the evolution of the plane-wave pertur- 
bations is sufficient for constructing the CMB transfer 
functions, summarized in Appendix B. 



For the plane- wave perturbations of intensity t Qj g(n), 
we use the standard basis vectors eg and e^, both or- 
thogonal to n, where eg lies in the (fe,n) plain and 
is orthogonal to it. In the basis {e.g,e^), a scalar (axi- 
ally symmetric about k) perturbation u a p has the general 
form [93] 

m-K'U'.-.)- < A18 > 

Its Stokes' polarization parameters are Q = q(x^,k-h) 
and V = U = 0. 

The photon polarization q 1 {x IJ ' 1 ni) evolves according 
to a transport equation [10, 14, 54] 

<? 7 + n l \I l q 1 = — {-q 7 + [(n.V,) 2 - V 2 ]tt p }. (A19) 

In the last terms of eqs. (A19) and (A13), 

7T p = C7 7 + p 7 o + Pj2- (A20) 

The variables p 7 o and p 7 2 arc defined by the following 
expansion of q^ijii) into scalar multipole potentials [69]: 

|g 7 K) = -V^o + SfaiV^V 2 ^! + (A21) 

oo 

+ ^(-l)<(2Z + l)P ( (n;f)VV 

1=2 

The hierarchy of evolution equations for the polarization 
multipole potentials follows from eq. (A19) as 15 

P-yO = V 2 p 7 i - i (p 7 o - \ 1tp) , 

~ 1 ~ 2 1 ~ 

Pyl = 3P7O — 3P72 — — P-yl, 

(A22) 

P 7 2 = -f V 2 p 7 l + § V 2 p 7j3 - i (p 7 2 - ^ 7T p ) , 



d. Baryons 

The evolution equations for baryons are similar to 
eq. (Al) for CDM. The growth of the baryon density 
perturbation is, again, minus the divergence of their ve- 
locity. The equation for baryon acceleration contains all 
the terms of eq. (Al) and an additional photon pressure 
term, related to the last term in it 7 equation (A15) by 



In momentum space, eqs. (A15) are related to the analogous 
Newtonian-gauge equations of Refs. [10] and [15] by elementary 

substitutions 8 1 = 4Aftf — ► 4(|<i 7 + *), 7 = 3kA^ -* k 2 Uj, 

and F 7i = 4A^p -> ffc'dyj. 



5 In momentum space, eqs. (A22) are related to the equations of 
Refs. [10] and [15] by G 7 ; = 4A^> -* |fc i+2 p 7 , for I = 0, 1 and 
G 7i = 4A£p -» |fc'p 7i for i > 2. 
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local momentum conservation in photon-baryon scatter- 
ing [10]: 



(A23) 



u b = -Hu b + $ - (Ub - Uj) . 
4. Quintessence 



Quintessence [30, 31, 61, 94] is defined as a classi- 
cal scalar field with canonical Lagrangian density C = 
K — V(4>), where K = — ^g^ v (f> 4l (f)^. In contrast to the 
previously considered species, the pressure of a classical 
field is not a local function of the field energy density. As 
a consequence, nonadiabatic pressure and the right-hand 
side of the density conservation equation (19) are now 
not zero. Nevertheless, if the field and the species which 
couple to it are perturbed internally adiabatically then, 
according to Sec. IV A, the field density perturbation d^ 
is time-independent on superhorizon scales. 

The evolution of quintessence can be described by the 
field equation (f>i^ = V' (</>), linearization of which in the 
Newtonian gauge yields 



2HS<f) + (k 2 + a 2 V")5(j) = 
= <f> ($ + 3*) - 2a 2 V'$. 



(A24) 



The Newtonian potentials $ and \& are affected by both 
the field perturbation 8(f) and its rate of change 8(f) 
through the field energy-momentum tensor Tj^} = 

When the scalar field is the dominant component of the 
universe, e.g., during inflation, the Newtonian potentials 
are equal and unambiguously determined by 8<p an d 8(f). 
Then 8(f) and $ = ^ can be elegantly combined into a 
Mukhanov-Sasaki variable v — a(8(f> + ^$), obeying a 
simple closed equation of an oscillator with a varying 
mass [9, 95, 96]. 

However, after reheating a scalar field no longer dom- 
inates the universe energy density. In this case the 
field dynamics can be integrated, at least numerically 
in the Newtonian gauge using the general equations (19) 
and (31). 

In the zeroth and linear orders of perturbation theory 
p 4 = K + V and P4> = K - V, with K = -±g 00 <f) 2 , and 
T-° = d G (j)di(j). Hence, the variables (d^, u$) are connected 
with (8<p, 8(f) as 



= -r-+a 2 V'^--$-3*, 
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Since, by the last equation, Sp,p = 8p$ — 2Vu^ and since 
anisotropic stress of a minimally coupled scalar field is 
zero, eqs. (19) and (31) yield: 



u,j> 



s(c) 



(A26) 



Here, 8^ = d$ + 3Hu^ + 3^ [it equals Sp^/ip^ + p^) 
in the gauge comoving to the field]. The ratio V/K = 
2a 2 V (4>)/<f> in eq. (A26) is set by a quintessence back- 
ground "equation of state" = p^j 



K-V . 
K+V ■ 



^ = -3W(l-c 2 ad ), 



J (j> ad 



- P<P — 



i; = w <t> 



3(l+iu^,)W 



Note that the "adiabatic sound speed" 



i fid 



does not 



describe any physical propagation. Inhomogencities of (f> 
propagate at the speed of light [1, 94], as evident directly 
from the field equation (A24). 



All 



In each of the above cases the perturbative dynam- 
ics in the fixed Newtonian gauge is formulated as an 
explicit initial-value Cauchy problem. For the consid- 
ered matter content, an arbitrary set of initial values 
for (d c ,u c ,d v i,dyt,P'yl,db,Ub,d^,u^) fully specifies phys- 
ically distinct and consistent initial conditions for scalar 
perturbations. 



APPENDIX B: CMB TRANSFER FUNCTIONS 
AND CVS 

Here we summarize the formulas for the angular spec- 
tra of CMB temperature anisotropy and polarization pro- 
duced by scalar perturbations. To be concrete, we con- 
sider adiabatic primordial perturbations with the primor- 
dial power spectrum 

(Cin(fe) &(*')) = (2nrS^(k-k')P c (k). (Bl) 

Correspondingly, we normalize the perturbation vari- 
ables in the transfer functions below [eq. (B4)] to Cin = 



CI 



fe/W(r)<l 



1. For the general nonadiabatic scenarios, 



P((k) would be replaced by a matrix of the primordial 
correlations of all the possible scalar modes with the mo- 
mentum k. 

The angular correlation spectra of the observed CMB 
temperature and polarization are 



*X Y \ 



2k 2 dk 



T*{k)T?{k)P c {k). (B2) 



In this equation the superscripts X and Y denote ei- 
ther temperature anisotropy (T) or curl-free polariza- 



(A25) tion (E) [16, 97, 98], and are the corresponding 

observed multipoles. The transfer functions T l x (k) are 
given by the line-of-sight integrals [15, 16] 



Tf> E {k) 



drS T ^ E {T,k) 3l {k{T -T)), (B3) 



where ji(x) = (—x) 1 (jgj) is a spherical Bessel func- 
tion, and the sources of temperature perturbation and 
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polarization equal 
S T (r,k) = g 



0_ 

dr 



©off + U b h - 7T p ( + — 



d 2 _ 

i 



l_ 2 

3 



111 cqs. (B4), ©off is defined by eq. (48), w p by eq. (A20), 
and 



g ($ + * 



(B4) 



fli(r) = exp - / 



(B5) 



S B (r,fc) 



+ 2)! 
0-2)! 



1/2 



4(r - r)2 



is the probability for a CMB photon emitted at time r 
to reach the Earth unscattered. 
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